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Abstract. We present new nonlinear differential equations for spacetime correlation 
functions of Fermi gas in one spatial dimension. The correlation functions we consider 
describe non-stationary processes out of equilibrium. The equations we obtain are 
integrable equations. They generalize known nonlinear differential equations for 
correlation functions at equilibrium [1] [2l IS [4] and provide vital tools to study non- 
equilibrium dynamics of electronic systems. The method we developed is based only 
on Wick's theorem and the hydrodynamic description of the Fermi gas. Differential 
equations appear directly in bilinear form. 
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1. Introduction 

Various spacetime correlation functions of free Fermi gas in one spatial dimension 
(and models related to it - impenetrable Bose and Fermi gases, XF-spin chain, 2D- 
Ising model, etc.) are Fredholm determinants or their minors and, therefore, generally 
speaking, obey integrable nonlinear differential equations with respect to spacetime and 
other parameters like temperature, chemical potential, etc. In practice, a derivation 
of these equations is a complicated task. Few equations are known. However, once 
obtained, they are indispensable in studies of off-shell properties of Fermi gas. 

Historically nonlinear differential equations for correlation functions first appeared 
in studies of the 2D Ising model [Ij. In Ref. [2] these equations has been extended, 
applied to the XY-spin chain, and have been derived in explicitly bilinear form with 
the help of Wick theorem generalized for these purposes in Ref. |5]. In Refs. [3] 
it has been showed that at zero temperature the two-point, equal-time correlation 
function of impenetrable bosons and of the XF-spin chain can be expressed through 
Painleve transcendents. These equations have been extended to time- and temperature- 
dependent correlation functions at equilibrium for impenetrable bosons and the XY- 
chain in Refs. [U [6] . 

The derivation of the known equations was essentially based on equilibrium and 
translational invariance properties of correlation functions. Although they have been 
recognized as integrable, their relations to integrable hierarchies remained unclear. 

Recently, interest focused on electronic systems out of equilibrium [7j and especially 
on non-stationary properties of propagating localized states [8]. There the physical 
properties change drastically and the dynamics acquires essentially nonlinear features 
such as hydrodynamic instabilities [9l [H]. Apart from being of a fundamental interest, 
research in non-equilibrium degenerate Fermi gas is driven by a quest to transmit 
quantum information in electronic nanodevices and to control and manipulate entangled 
quantum many-body states. Integrable equations describing the non-equilibrium - non- 
stationary states of electronic gas, once available, would be a valuable tool in studies of 
a complex dynamics of electronic gas. They are derived in this letter. 

We will obtain nonlinear differential equations with respect to spacetime for the 
correlation functions of vertex operators 

where ^Pr^l are right/left handed chiral Bose fields of fermionic currents and aL^pi are 
arbitrary (not necessarily real) parameters. The definition of these fields in terms of free 
fermions is given below. Equations with respect to parameters of the density matrix 
(e.g., temperature or chemical potential) can be also obtained, but are not discussed 
here. 

The correlation functions of vertex operators ([1]) are important in many 
applications. In the case of equilibrium they also satisfy some previously known 
equations. A few particular cases are worth mentioning. The case ai = = 1/2 
describes a one-point function of impenetrable bosons [31 H] . The imaginary ai = 
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and = —af> yield a generating function of moments of the distribution function of 
the density and velocity of electronic gas. The chiral correlation functions a^, = or 
aji = describe a tunneling into an electron gas (Fermi Edge Singularity) ^0\, etc. 

In this article we give a simple physical derivation of equations on these correlations 
functions. Our derivation is based solely on a general form of the Wick theorem 
[llj , the representation of Fermi-operators through Bose fields (see e.g., |12j), and 
the hydrodynamic description of Fermi gas [13j. Our equations appear directly in 
Hirota's bilinear form [14j. They are known as modified Kadomtsev-Petviashvili (or 
mKP) equations and are the first and the second equation of mKP hierarchy. In a 
particular case, in free fermion systems these equations have previously appeared in 
[T5| |9]. The nonlinear equations, being integrable are proven to be an effective tool in 
computing asymptotes of correlation functions in different regimes. These computations 
are specific to each particular problem and will not be discussed here (see, Refs. P E\ 
for some practical applications). 

In the following after necessary preliminaries on the Fermi-Bose correspondence, we 
describe the correlation functions ([5]), present the nonlinear differential equations they 
obey ( fT6]ll8|ll9l) . and then give a short derivation of those equations. 




where ipp is a mode of a fermion field ipi^x) = J^p^^^^i^p- 

A central point of our approach is the hydrodynamic description of the Fermi gas. 
We briefly review it. The properties of the Fermi gas are fully described in terms of 
canonical hydrodynamic variables: density and velocity 




2. Fermions and a Bose Field 



We consider free fermions on a circle of a circumference L 



(2) 




where the chiral (right, left) currents are 



iR,L/ \ '2,11 ^ m ^ X 




(3) 



lead to a canonical relation between density and velocity: 

h 

[p{x,t), v{y,t)] = -i—V6{x-y) 



(4) 



and to the canonical Bose field 



VipR,L = ±J^'^(a;, t) = j-v± np, [cpL, (pn] = 0. 
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We will be interested in the spacetime dependence of correlations of two vertex operators: 

Here g is the density matrix and ^ = (x,t). The colon denotes normal ordering with 
respect to the vacuum [16]. 

Below we drop the chiral subscripts R, L in all formulas and assume that the upper 
(lower) sign corresponds to the right (left) sectors. 



3. Coherent states 



We will consider correlation functions with respect to coherent states. This means that 
the density matrix is an element of gl{oo), i.e., it is an exponent bilinear in fermionic 
modes: 

Vp.? / 

This choice of the density matrix allows us to use Wick's theorem. 

The set of coherent states is rather general. It exhausts most interesting 
applications. For example, the density matrix can be a Boltzmann distribution, or any 
other non-equilibrium distribution as in f7]. Another choice would be ^ = where 
\V) = exp (J2k tkJk) |0) is a coherent state (i.e., an eigenstate of current operators). Such 
states appear as a result of shake up of the electronic gas by electromagnetic potential 
with harmonics t^. 



4. Quantum Hydrodynamics 

The Hamiltonian of free fermions ([2]) in the sector of coherent states can be expressed 
solely in terms of density and velocity [13] 

^ rf^^3A,,. (6) 



2 6m 
Evolution of density and velocity follows 



: p : +V : (pv) := Continuity equation 



: V : +- V : (v^ H ^P^) := Euler's equation 

2 

The hydrodynamic equations are decoupled into two independent chiral quantum 
Riemann equations 
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5. Vertex operator and correlation functions 

The Heisenberg evolution equations for the chiral vertex operators follow from the 



quantum Riemann equation ([7]) (see [Appendix A\j 



(idt- —V''):e'^^: = —a(a + l):e'^'^T:, (8) 

+ — V'):e*("+^)^: = - — a(a + 1) :e*('^+^)^f :, (9) 
2m 2m 

where T = : (Vy?)^ : —iV'^ip and T = : (Vy?)^ : +iV'^ip are holomorphic (antiholomorphic) 

components of the stress-energy tensor of a chiral Bose field (with the central charge 

1/2). 

A particular consequence of these equations is that the vertex operators e^**^-^'^ 
obey the Schroedinger equation. This yields to a familiar formula representing fermions 
as exponents of a Bose field 

-.e'"^:- ^/L^p{x), : e-'^^ -.^ ^i)\x) . (10) 

A useful representation for current J = ±: e"*'^(— iV)e*''' : and the stress energy tensor 
reads 

T = - : e^^'^V^e^'^ :, f = - : e^'^V^e"*'^ : (11) 

Current and stress energy tensor can be further cast in fermionic form using ( fTOj) and 
OPE for vertex operators as follows from the current algebra pil|T7] 



With the help of ([T2l [TTj) and ([TO]) we write 

T = ^4mV"^VV', f = ^47^#VV'^ J = 27r^V- (13) 

Finally we cast the equation of motion of the vertex operator in a suggestive bilinear 
form (Quantum Hirota equation) 

( tDt - —Dl) : e"^ ■ 6'^"+^^^:= 0. (14) 
Y 2m J 

Here D is the Hirota derivative Df ■ g = dfg — fdg and D^f ■ g = d^fg — 2dfdg + fd'^g. 
This equation follows from the quantum Riemann equation ([7]). 

6. Nonlinear differential equations 

The nonlinear equations themselves do not depend on the parameters aL,R, hL,R in ([5]). 
For brevity we put L = 27r in formulas below and consider only the chiral and neutral 
(a = h) correlation function, which we denote as 

r,(e,0 = tr(:e-^'^^«')::e^'^^«):f?). (15) 
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Being presented in bilinear form the differential equations for correlation functions 
look identical to the quantum equation ([Ti 




'a+l 



0. (16) 



The equations must be viewed as equations for ra+„(a;, t; x', t') in the continuous variables 
X, t and the discrete variable n, where \a\ < 1/2. As such they are a closed set of 
equations known as the mKP equations. The correlation function can be found given 
initial conditions at t = 0, for all integer n and any x. Certain asymptotes or analytical 
conditions may also be helpful in solving the problem for particular cases. This program 
is taken up in practice in [9l [8], where in [9] the density matrix had the form |v4)(0| 
where |0) is the ground state of electron gas and \A) is a chiral coherent state of the 
form \A) = g\0). We also comment that where \a\ < 1/2 is a Fredholm determinant, 
while the Ta+n are its minors. 

The mKP equations (fT6l) has many different forms discussed in the literature. We 
will mention one of them. Setting 

$ = log , $ = — log TaT^+i 

m \Ta+i) m 

the equation f|T6l) becomes 

9i<l> = -^(V$)2 + Av^l.. (17) 
2 2m 

In this form the classical equation (ITTl) resembles the quantum equation ([7]). The last 
term in (1171) can be considered as a quantum correction to the semiclassical Riemann 
equations. 

In addition to the evolution f|T6|) there is another equation which does not involves 
time evolution. This equation connects correlation functions with different density 
matrices: g and g = Y.pqO'pq'Ppff'Pqi where apq is an arbitrary non-degenerate matrix. 
Then and fa = tr e'^'^'^^?') ::e'''^(«) : ^) are related by 

DxD^'fa ■ Ta = a^{fa+lTa-l + r„+if„_i). (18) 

This equation is a minor generalization of familiar 2D Toda lattice equation. 

An important particular case occurs when the modes p,q in g lie very far from the 
Fermi surface. Then ~ and the equation flTSl) becomes ID Toda chain equation 

D^D^>Ta ■ Ta = 2a^Ta+lTa-l. (19) 

This equation, together with f|T6|) being written for pairs Ta,Ta+i and Ta-i,Ta give a 
closed set which involves only ra_i, r^, Ta+i- 

7. Translation invariant states 

Further closure can be achieved for special cases. For example, if the density matrix 
commutes with momentum, the correlation function depends on the difference x — x' 
(e.g., the density matrix depends only on occupation numbers rip = iphpp, g = 
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exp^^pApTip. Then it also commutes with the Hamiltonian, is a function of t — t', 
and describes stationary processes). In this case we denote the tau-function as 

'^a (-^7 X ^t ^ Fa (-^ X ^ t ^ t f 

As follows from (I16fl9p three functions Fa-i, Fa, Fa+i obey three equations with respect 
to (x, t) 



iDt + 1 Fa ■ Fa-l = 0, 



DlFa ■ Fa = -2a^Fa-iFa+i. (20) 

Here we used D^' = —D^. 

These equations are the bilinear form of the Nonlinear Schroedinger equation 
(without a complex involution). Introducing \E' = Fa+i/Fa and ^ = Fa-i/Fa we have 

\ 2m J m 

y 2m J m 

These equations were obtained for time and temperature dependent two-point 
correlation function for impenetrable bosons at equilibrium (g = e~^^'^) [4j. The 
impenetrable bosons are equivalent to free fermions. The creation operator of the boson 
is the vertex operator ([1]) with = = 1/2. 



8. Derivation of Nonlinear equations 

The fact that the quantum equation for vertex operators 0141] . and the classical equations 
for the correlation functions, (fT6l) . have the same form may not be an accident but may 
reflect a general property of integrable equations written in bilinear form. Here we limit 
our discussion by demonstrating this phenomena for the equation (ITGll . 
We prove now the main formula f|T6l) . First we adopt the notation: 

{Oi{x),02{x')) = {Ta{x,t;x',t'))-' 

X tr [g e'^' : e'"'^(")Ci(x) : e-'^^*"*'); e-'"'^("')C2(x') : e"'^*') , 

where only the spatial dependence has been specified in the r.h.s.. 

Let us muhiply ([8]) and ([9]) by : e'^^^^^'''*') : and : e-*('^+i)^(^''*') : respectively, take 
the trace, and substitute into Eq. flTB]) . We find that flTB]) holds due to the identity 

(e-^^, e'^f) + (1, T) (e-^^, e'^) = 2(1, J) (e"*'^, e'^ J), (22) 

where 1 is the identity operator, and we further suppressed the arguments x,x'. 

To prove this identity we first write : e*'^T : in terms of fermions. Using (iTTl) we 
write : e*'^T := — lim^^^^a, : e*'^*^^^e*'^*^^^e~*'^*^^) :, where we split spatial points. We now 
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use the formula (1121) for vertex operators and Eq. (ITT]) to write 

:e*^(")T(x): = Tii'^'nf'^ lim v 2 (l/ ~ "^^^^ ~ ''^ (xUlzUUy). (23) 



This formula allows to write the first term in (1221) in terms of a correlator with four 
fermions insertion 

T i{27^f lim V; ~ 7 (^t(^/)^ ^(3,)^(^)^t(^)). (24) 

Now we use a general form of the Wick theorem [TT] to express the four fermion correlator 
in terms of correlators containing two-fermions 

(^t(a;'),^(a;)^(;^)^t(^)) = (^(a;),^t(^'))(^t(^)^(^)^l) 

- {ij{z),ij\x')){ij\y)ij{x),l). (25) 

After this, we trace back the path which lead us to (12^ . With the help of Eqs. (ITUl) 
and (fT2|) we obtain the formulas 

(l,V^(z)^t(^)) = ± __I__(i,e^^We-^(^)). (26) 

We use them in order to write the r.h.s. of (125!) in terms of bosons. After that, one 
applies the operation lim^^^_>.;. fa~^K!/-^) ^j^g (l25l) . takes the derivative, and 



merges the points. This leads to (|22|) . The calculations are somewhat cumbersome, but 
straightforward [18| . 

We briefly comment on the proof of (119!) . the proof of the more general (fT8l) goes 
along the same lines. First we write (1191) in the form 



(J, J) - (J,l)(l, J) = (e^^e-^^)(e-^^e'^). (27) 

The proof of this equation is easier: after replacing bosonic exponentials and the current 
operator by fermionic operators according to (I10|13l) one recognizes Wick's theorem. 



9. Integrable hierarchy of nonlinear equations 

As a final remark we mention that higher equations of the mKP hierarchy describe 
evolution of Fermi gas with an arbitrary spectrum H = J^p^pi^li^p- They are generated 
by the Hirota equation described in [15]. Equations for magnetic chains equivalent to 
the Fermi gas can be obtained in a similar manner. 
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Appendix A. Equation of motion for vertex operators 



In this Appendix we sketch the derivation of equations ( l8|9l) from the Riemann equation 
([7]). First of all we show how to calculate the time derivative of the vertex operator. A 
subtlety here is that ip does not commute with ip. 

We recall a simple consequence of the Hadamard lemma. If H and A are two 
operators such that [[H, A] , A] commutes with A, then 



H, e' 



E,A\^-{{R,A\,A\ 



\H,A\-)^^\\E,A\,A\]e^ 



If is a Hamiltonian dtA ~ [if, A\ and we obtain for the time derivative of an exponent 



2! 



A. A 



A 



2! 



A. A 



(A.l) 



Let us now compute an evolution of e*""^. We assume that <y9 is a right chiral field. The 
left chiral field obeys the identical equation. 

We notice that the commutator , vanishes, and then apply (1A.1I) . We 

obtain 



a : e*"'^ : 



2 



+ e 



ia<p ' 



iacp + 



laif — 



[la) 

"2r 

(ia) 
~2[ 



laip 



Let us now compute the commutator Using ([7]) and 

Fourier components 



m 



(A.2) 
we obtain for 

(A.3) 



Then we proceed as 

0+(a;),v3"*"(x) 



1 C^'^ 

m\L 

1 /2tx 
m\L 



p,q<0 



fc<0 fe<p<0 



1 27r 
hm L 



fc<0 



h 



m 



Repeating the same calculation for (f we obtain 

h 



m 



(A.4) 



Finally, using ([0]) in f[OD we obtain (TO 

[iaY h 



a : e*"'^ : 



h 



2 m 



^a^p 



taip 



(A.5) 
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